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Abstract. We study the problem of forecasting volatility for the multifractal 
random walk model. In order to avoid the ill posed problem of estimating the 
correlation length T of the model, we introduce a limiting object defined in a 
quotient space; formally, this object is an infinite range logvolatility. For this object 
and the non limiting object, we obtain precise prediction formulas and we apply 
them to the problem of forecasting volatility and pricing options with the MRW 
model in the absence of a reliable estimate of a and T. 
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1. Introduction 

In recent years, the Multifractal Random Walk (MRW) model introduced by 
Bacry, Delour and Muzy ([1]) has received much attention from the financial practi- 
tioners. The MRW appears as a natural extension of the basic geometric Brownian 
model (GB). Let St be the price of an asset; in the GB model we write: 

St = Soe^\ (1.1) 

where Yt = (yU — |o"^)t -|- crBt. Here fi is the mean return rate , a the volatility 
of the asset supposed to be constant, Bt a standard Brownian motion. The GB 
model plays a fundamental role in finance as it enables to give a price to options (by 
the famous Black-Scholes formula). Nevertheless when compared with reality the 
model displays severe drawbacks; more precizely, the GB model does not display 
the following stylized facts which are largely acknowledged in the litterature ([7], 

• The volatility fluctuates randomly and follows approximately a lognormal 
distribution. 
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• While the returns are rapidly decorrelated, the volatility exhibits long range 
correlations following a power law 

• The returns are heavy tailed. 

The above stylized facts lead naturally to propose as model of volatility a random 
measure called the limit lognormal model. This object, introduced by Mandelbrot 
([11]) in the context of turbulence, was rigorously defined and studied by Kahane in 
[To] under the name of gaussian multiplicative chaos. More precisely, Kahane devel- 
opped a general theory for cx-positive kernels (see [10] for the exact definition) and 
suggested how to use this framework to define the object introduced by Mandelbrot 
in [n]. 

Let us give a very sketchy presentation of the gaussian multiplicative chaos. Let 
T and 7 be given parameters (T is the correlation length). We consider the positive 
kernel p{s,t) = ln'''(T/|t — s\). To this kernel, we associate the stationary gaussian 
process Xj' with covariance p{s,t) (in a generalized sense since p{t,t) = 00). The 
associated multiplicative chaos is a random measure defined (formally) by: 

m{dt) = e^^"-^^«-'^")')rft. 

Of course, this formula has no meaning since E{{X'[Y) = 00 but it clearly suggests 
the limit procedure by which we can define rigorously the random measure (in [3], 
it is proven that the kernel p is cr-positive in the sense of Kahane and therefore the 
measure m is a particular case of the general theory developped in [lOj). 

Following an idea of Mandelbrot and Taylor ([I2]), one can consider the model (for 
the log price) of a Brownian motion subordinated by an independent random mesure. 
This was first performed with the measure m by the authors of [3]: this defines the 
MRW model. More precisely, we define the random time change d{t) = m([0,t]). 
The MRW model is then the four parameter stochastic process given by: 

Y^^ = pt-^a'e{t) + crBeit), 

where the Brownian motion B is independent of m. If t is not too large (a few years 
for example), the drift term pt — ^a'^6{t) is in practice negligible compared to the 
Brownian term so the MRW reduces to the three parameter [a, 7, T) process: 

From classical properties of the chaos m, it follows that 6{t) is a continuous process 
so that Y^'" is also continuous. By using the scale invariance of B, one can write (at 
least formally) as a stochastic integral: 

i;^ = a re^^"-^'^«^")')d5„ (1.2) 
Jo 

where A = | is the intermittency coefficient. One can then obtain Y^'" rigorously 
as the limit of a discretized approximation scheme; more precisely, let r be some 
positive step parameter, {en)nez be a standard gaussian i.i.d. sequence and (X^)„gz 
be a zero mean stationary gaussian sequence independant of e with kernel (see lemma 
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( 16. ip of the appendix for the existence of X^): 

If 4A^ < 1, one can show the convergence in law (in a functional sense) as r goes to 
of (Fi^'^)t ^0 given by: 

Vt/r\ 
n=l 

towards iY^)t^Q iM)- 

To what extent is relevant to model financial markets? As it gives a fat tailed 
distribution for the returns, long range correlations and a lognormal fluctuating 
volatility, the MRW appears as a rather relevant model for exchange markets. The 
main limitation of this model is the symmetry of the distribution of returns and 
therefore it can not capture the leverage effect observed on stocks and indices (see 
[5] for a quantitative study of the leverage effect) . One of the most important 
predictions of the process {Y^)t^o concerns the logvariogram. More specifically, 

let Y'^ = max Y,^ — min Y,'^ . Then, one defines the logvariogram V(j) by the 

^ te[j;j+i] te[i;j+i] 

formula: 

vu) = E[{HY;')-HY,^)n 

One can show that there exists C > such that for all j ^ T: 

VU) = C + 2X'Hj) + X'6U,X') (1.4) 
where C ~ 0.29 and the terms 6{j, A^) do not depend on T or cr such that: 

sup 5(j, A^) -> 0. 

In practice, one can neglect the last term of the above expansion since the value of 
A^ for a financial asset (stock, currency, index, etc..) will typically belong to the 
interval [0.01,0.06]. 

We see on the figures [1] and [2] the empirical logvariogram of the SP500 index 
on three disjoint periods of 6 years and the empirical logvariogram of the currency 
Euro/ Australian dollar on the period 2001 — 2007 along with the corresponding 
regressions of the form C + 2A^ln(j). One can notice that the logvariograms of the 
SP500 are less and less noisy as time evolves (perhaps due to the increase of liquidity: 
tick size, volume, etc.). We propose to estimate the intermittency parameter A^ by 
the aformentionned regression: Monte Carlo simulations performed with 1000 trials 
of a MRW on a period of 6 years with A^ = 0.02 show that the above estimator 
belongs to the interval [0.01,0.03] with a confidence interval of 90 percent. 

Finally, we mention the problem of the estimation of the parameter T. Typically, 
T is very hard to estimate precisely (With the financial data available, one finds 
huge error bars on the estimation of this parameter: for typical estimates of T, see 
[2]). Roughly, T will exceed 2-3 years for an asset. It is easy to convince oneself 
that beyond a window of size 5 — 6 years, one can not assume that the sequence of 
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returns of an asset forms a stationary sequence because markets evolve in time. In 
these conditions, ergodicity of the data breaks down and one can not rely on a precise 
estimate of T. In fact, it is even impossible to determine T with data on a window 
of size T; indeed, if Y'^ and Y'^ are two MRW with respective parameters (a. A, T) 
and {a, A, T) such that T ^ T then the following generalized scale invariance holds: 



with parameters (^ae^^~^^^''^^'^^^\ A,T). In conclusion, the statistical estimation of T 
and a is an ill posed problem and the most natural way to get rid of the problem 
is thus to let T go to infinity and work with the limiting object. In doing such a 
procedure, one must work in a quotient space (see definition (12. 3p where one views 
the log volatility X as an element of the quotient space i5'(M)/R). 

The rest of the article is organized as follows: in section 2, we give some preliminary 
definitions: in particular, we introduce the log volatility processes and X. In 
section 3, we give explicit prediction formulas on X"^ and X. In section 4, we show 
how to use the results of section 3 to forecast volatility and price options with the 
MRW model with no knowledge on the value of a or T. Finally, in section 5, we give 
the detailed proof of the results of section 3. 



2.1. Definitions. For a < 6, we denote by V{]a,b[) the space of smooth functions 
on ]a,b[ with compact support included in ]a,b[ and by V'(]a,b[) the space of dis- 
tributions on ]a, b[. Let S{R) denote the space of Schwartz functions and iS'(M) it's 
topological dual (the space of tempered distributions). The Fourier transform of an 
element / of iS'(M) is defined as usual by: 



We remind the definition of the classical Sobolev Spaces H^{M.) for s G M: 




2. Definitions and preliminaries 



V0e5(M), </,0>=</,0>, 



where is the classical Fourier transform of 0: 




H%R) = {/ G S\R) ; feL 



1 

7oc5 




The space H^{R) equipped with the scalar product: 



<f,9>H^= [ (i+ier)7(oim 
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We also introduce the homogeneous Sobolev Spaces H^iM.) for s < i; 

W{R) = {/ G S'{R) ; / G LL, / l^nhOl'd^ < oo}, 

Jr 

The space 7i*(M) equipped with the scalar product: 

Jr 

is a Hilbert space. 

For s G [|; |[, the homogeneous Sobolev Spaces H^iM) are defined by: 

n'{R) = {f e S'{R)/R ; f G n'^\R)}. 
For s G [|; |[, the space 7i''(M) equipped with the scalar product: 

1 



is a Hilbert space. 



2.2. The processes (Xj)t6M, (Xt)fgM and their reproducing kernel Hilbert 

space. We consider the reader is familiar with the theory of gaussian processes. For 
a complete account on the theory, we refer to [1]. 

We introduce the generalized gaussian process {Xj)t(z^ as a gaussian measure on 
the space iS'(M) with the following characteristic functional: 

V0 G 5(M), E{e'^^^^'^^^'^') = e-^/m(k4>(t)xTdt)^) 

= e-V2/ /b2 <t>m{s) \n+{T/\t-s\)dsdt_ 

It follows from Minlos's theorem that the above formula defines a unique gaussian 
mesure on iS'(M) ([9], chapter 4). Let sine be defined by the classical formula: 

. , , , sin(t) 
smc X = / at. 



t 

By using Parseval's identity and lemma EH] of the appendix, one can deduce that: 

1 /■ sinc(27r|^|T)| -j^^ 2 




0(t)0(s) ln+(T/|t - s\)dsdt = - / V \mm- (2.1) 

Jr KI 

The reproducing kernel Hilbert space Ht ([!]) associated to {Xj')t^^ is therefore 
the space: 

HT = {feS'm; [ ■ rJ^' |^J /(OP^e<oo}. 

Jr smc(27r|^|T) 

Since there exists two constants c, C > such that: 

the space Ht is the space if^/^(M) equipped with an equivalent norm. 
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One would like to let T go to infinity in the above definition. In order to give a 
rigorous meaning to that procedure, one must work in the space iSo(M) of Schwartz 
functions of average zero: 



5o(M) = {0 G S{R) ; / (f){t)dt = 0}. 

Jr 




The topological dual of iSo(M) is then the quotient space 5'(M)/]R. We therefore 
introduce the generalized gaussian process {Xt)t^^ as a gaussian measure on the 
space iS'(]R)/M with the following characteristic function: 

= e-i/2//B.2 0(t)</>(s)in(i/|t-s|)dsdf_ 1^2.4) 
By letting T go to infinity in identity (12. ip . one gets: 

<pm{s)Hi/\t-s\)dsdt = l [ ^ll^rfe 

Jr \^\ 

The reproducing kernel Hilbert space H associated to (XjjgK is therefore the space: 

H = {fe S'{R)/R ; 2 / 1^11/(01'^^ < oo}. 

Jr 

The space H is thus precisely the space 7Y^/^(M). 

Notations: In the sequel, we will use the notations kx = ln'^(T/|.|) and k = 
ln(l/|.|). 

3. The prediction formulas for X,X'^ 

3.1. Prediction formulas for (XjfgK. We have the following explicit formula for 
the conditional expectation of {Xt)t^^: 

Theorem 3.1. Let L be some finite positive real number and fit) the trace on 
the interval ] — 2L,0[ of a function in 1-(}^'^{R^. Then, the conditional expectation 
E[Xt\{Xs)^2L<s<o = /] is a function ofH^^'^iR) that satisfies: 

FIXKX) -fl- ^ r ^-9L{tY f{s)ds 

T^L J_2L (1 - 9L{t)y - 2gL{t)j^ ^1 - (1 + f)2 

where the application gi is given by: 



gLit) = ^ + l-Jij + ir-l, zft>0, 



or 



gLit) = ^ + l + ^{^ + iy-l, zft<-2L. 

Remark 3.2. Since a function in 'H^^'^(K) is locally in for all q < oo, we easily 
check that the integral in \3.1\ makes sense. 

By letting L ^ oo in the above formula, we get the following corollary: 
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Vt>0 E[Xt\{Xs)s<o = f] = - I ^ ^ f{s)ds. (3.2) 



Corollary 3.3. Let f{t) be the trace on the interval ] — oo,0[ of a function in 
7Y^/^(M). Then, the conditional expectation E[Xt\{Xs)s^Q = f] is a function of 
T-C^/^Ir) that satisfies: 

The proofs of theorem I3.H corollary 13.31 and theorem 13.51 (see below) are given 
in section 5. In the sequel, we will note KL{t, s) the kernel that appears in formula 
fl3.ip and K{t, s) the kernel that appears in formula 03.21) . 

Remark 3.4. Since (Xt)tm does not belong almost surely to 7i^/^(M), we remind in 
what sense one must understand formula liS. 1\) (the same remark applies to formula 
( 15*. ^) ). An equivalent formulation of 1\) is for all in Sq(R) with support in 
] - cx),-2L[U]0,oo[ : 

E[ [ <p{t)Xtdt I ( [ tlj{s)XJs] ] = [ ( [ (P{t)KL{t,s)dt] XJs 

\J~2L / ^eV{]~2Lfi[);f ip=0 J -"iL \Jr J 

3.2. Consequence on the prediction of (X'^\^^. As a consequence of the above 
formulas on [X^i^^, it is possible to obtain an exact prediction formula on (X^^)tgK.This 
formula is no longer defined modulo constants and is valid whatever the 
value of the correlation length T . It is obtained as a perturbation of the case 
T = oo. 

In order to state the formula, we introduce Lp as the unique distribution of iir~^/^(M) 
with support in [—1,0] solution of: 

VtG]0, 1[, y ln(l/|t- s|)(^(s)ds = 1. 

We refer to section 5.3 for the existence and uniqueness of In fact, the distribution 
(y9 is a function and it is given explicitly by the following formula: 



V^'-^l-l.OI. ^W = ^^=i= (3^3) 



This leads to: 



Theorem 3.5. Let L be some finite and positive real number such that 2L < T and 
f{t) the trace on]-2L,0[ of a function in H^/'^{R). Then E[X^\{XJ)^2L<s<o = f] 
is a function o/i7^/^(M) that satisfies: 

VtG]0,T-2L[, E[X^\{Xj)^,L<s<o = f]= r KlA^^ s)f{s)ds, (3.4) 
where the kernel satisfies for t in ]0,T — 2L[: 

A remark similar to remark 13.41 applies to formula (13.41) . 
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4. Application to volatility forecasting and option pricing 

In this section, we suppose that we ignore the values of a and T but 
that the value of is known. 

4.1. Forecasting volatility. 

Discretizing formula Ii3.2\) . Let us first suppose that T is large compared to L so 
that the prediction kernel is given by K^. Suppose also that the time t at which we 
intend to predict is small compared to L. It follows that we may use the prediction 
formula (lOl : 

E[X,|(X,),<o = /] = - r ^! ^ f{s)ds. 



71" J~oo {t - s)^/^' 

Let us take some small time lag r and discretize the above formula. We get: 

N 

/ ~ ^ l[.^k+l)r,-kr[f{-kT) (4.1) 
fc=0 

and thus for n ^ 1: 

^ N n—kr I 

E[X„^|(X_fc^)o^fc<;7V = (/(-A;r))o^fc^7v] ^ -V / 7 -^dsfi-kr) 



k=0 ' 

Let us denote 



J(k-l)/n (1 + S)^/S 



2., k , k-1, 

= — (Arctan-i / Arctan-i / ), 

TT \ n \ n 

so that the discretized formula writes: 

N 

E[XnT\{X_kr)o k ^ n] — y^Qn.fc+l^-fcr- 

k=0 



Of course, we have for all n ^ 1, Efc=o = ^^^^ctany'- so Efc=o<,fc+i - 1 

for n <^ A. Notice also that a* does not depend on r. 

Application of the discretized formula. We will model the (log) return (r„)„g2 at 
scale r of a financial asset discretized MRW: 

where (e„)„g2 is the noise process (i.i.d. standard gaussian) and {crn)nez the volatility 
process: 

with given by definition fll.3p . 

Suppose one can observe the historical volatility (T„ (or equivalently X^) on some 
time window {— A^, ...,0} (Of course this is rigorously not the case in finance and 
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one must use filtering theory, intraday data, etc... to get a proxy of (T„). We can 
decompose tfie process X^, n ^ 1: 



N 



= ^n,k+lX^k + (*) 



k=0 

wliere Z„ is a centered gaussian process independent from XI^,...,Xq and the 
coefficients an,k+i are uniquely determined by the conditions: 

i?[x;z„] = o, j = -N,...,o. 

If we suppose that N <^T/t and n <^ N then formula (*) appears as a discretization 
of formula (13. 2p and we may approximate: 

SO that ^fcLo'^n,fc+i — 1- Then we write: 

N 



fc=0 
N 



k=0 

from which we get the conditional expectation: 

N 
k=0 

Now let us calculate E[Zl]. We use: 

TV 
k=0 

and 

TV 



fc=0 

Therefore, we get: 

fOO POD 



1 fnr 

Y (y*n,k^xX-k - - / 7 y ^sds 

-n J-oo[nT - s)^s 



-1 /"OO noo -1 

= ln(n) + - / / H\s - ^)dsd~s 

TT^ Jo Jo (1 + s)(l + s)JsVs 



(1 + s){l + s)y/sVs 

If we set C equal to the constant term in the above expression, we get the following 
simple formula: 

TV 

E[an\{a.,)o^k^N]^e'^n''/'l[{a.k)<'^^K 

k=0 
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One can notice the remarkable fact that, in the above formula, ex, r, T have disap- 
peared. Nevertheless, the formula depends on and we get the following sensitivity 
with respect to A^: 

^^^^^^^^^^^i^ = (f + i^)BKI(.-.)„.««]. (4.2) 

Since C ^ 1.33 numerically, if one considers the problem of forecasting volatility 
on a period of 1 month (n ^ 20) and A^ is known with a precision 6X^ = 0.01, the 
above formula fl4.2p implies that the forecast is correct with a very high precision of 
approximately 2 percent. 

4.2. Option pricing with unknown parameters. Let St be the price process 
of a financial asset. We model the log price by the continuous MRW Y'^ (with 
underlying parameters (a. A, T)) given by expression (11.21) . In this subsection, we fix 
a real positive number t. Suppose one can observe the historical (log) volatility Xj 
on some time window ] — 2L, 0[. The price of a standard call option with maturity 
t and exercise price K is then given by (for simplicity, we suppose the risk free rate 
r = and that there are no dividends): 

E[(5oe^" - KU)\{X^).2L<s<o = /], (4.3) 

where Y^'" is given by the approximation (14.71) . For simplicity, we will also suppose 
that the maturity t is not very large (typically, a few months) such that one can 
make the following approximation by an additive model: 

5oe^" ^5o(l + Fi^). 

Then, one can perform a cumulant expansion in formula (14. 3 p up to the kurtosis (see 
p. 244-245 in ^ for details) and obtain the following smile formula for the implied 
volatility T,[K,t) as a function of the strike K and the maturity t: 

m.t) = %{l + nt{ ^^~^f -I)). (4.4) 

where is the forecasted variance: 

^l = E[{Yrf\{Xj)^2L<s<o = f] (4.5) 
and Kt the forecasted kurtosis: 

_ E{{Y;'Y\{XJ)^2L<s<q = f] _ ,^ 



We outline the computation of (14.50 . Let r be some small observation scale and let 

■r(n-l) 



— YJ',_^-. be the (log) return at scale r. We have: 



n=l 
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One then uses the expression r„ = ay/re^^" ''^*^^/'^)e„ of the previous subsection 
and the discretization: 

N 

f -Yl l[-(fc+l)T,-fcr[/(-fcr). 
k=0 

We therefore get with = [2L/rJ: 

n=l 

The computation of each term 

is similar to the computations of the previous subsection. 

Since, in practice, is roughly found around the value 0.02, we perform an 
expansion of for ^ 1 in order to compute Hf We have the following scaling 
identity: 

(Law) 

where e is a standard gaussian variable independent of 0{t). We can derive formally 
the following series of approximations for A^ -C 1 (we will use for <C 1 and X 
a generalized centered gaussian variable the approximation e'^^~'^ ~ 1 + cX: 

see [2] for an exact mathematical formulation): 







{l + 2\Xj)dse 
= aVt^ ^ + T / ^'^''^ ^ 

^av^eT/o^="'^-^^W^)e (4.7) 

where e is a standard gaussian variable independent of X'^. Using expression (14.71) . 
standard computations give the following expression for k^: 



If we let T go to infinity, we get the following limit expression: 

T^oo 
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Remark 4.1. The formulas ( [^.<^| j and ( [^.^^| ) are independant of f , which is a con- 
sequence of the lognormal assumption: every lognormal model for the volatility will 
give a constant conditional kurtosis. Intuitively, this can appear as a limitation of 
the model. 

Remark 4.2. It is remarkable that the conditional kurtosis given by formula Iji4-^ 
tends to a finite limit given by Iji4-^ T goes to infinity. Indeed, the unconditional 
kurtosis tends to infinity as T goes to infinity. 

5. Proof of the prediction formulas 

5.1. Preliminary results and definition. In this section, we prove the theorems 
of section [31 Let 11 denote the upper half plane Rx]0; +oo[ and V^^(n) denote the 
following quotient space: 

V\lVj = {ueV'iU)/R; [ \\/ u\'^ dxdy < oo} 

Ju 

equipped with the norm 




u\v^{u) = W / \'Vu\'^dxdy. 



n 



On the Hilbert space V^^(n), one can define in a classical way a trace operator: 

Lemma 5.1. There exists a unique continuous linear operator T : V\U) H^^^iR) 
that satisfies: 

v/g5(m2), r(4>o)(a^) = /(x,o). 

The operator T is onto and satisfies: 

|r(/)l^i/2(K) = ^1/1^1(0)- 

We introduce the Poisson kernel P by the following formula: 

V/ G n'/\R), P{f){x,y) = Jl'} f{z)dz. (5.1) 

Notice that the existence of the above integral is not obvious: one first defines it for 
/ in S(R) and then extends it to H^/'^iM) by density 
The next lemma relates P to a variational problem. 

Lemma 5.2. For all f E 'H^^'^{'R) there exists a unique solution to the variational 
problem: 



inf ( / \Vurdxdy). (5.2) 
«eyi(n);r(n)=/ Jn 

and the solution is P{f) given by the Poisson formula \5.1[ Moreover, P{f) is the 
unique solution to the Dirichlet problem: 



(5.3) 
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and the following identity holds: 

l^(/)lyHn) = 2vr|/liv.(M)- (5-4) 

Proof. Let / G 7Y^/^(M) be fixed. A simple computation shows that the M? Fourier 
transform of P{f) is given by: 

W)fe,6) = -/te)JT72- 

SI + S2 

Therefore one gets: 

\Pif)\vHu) = l [ \VPif)\'dxdy 



= 2n / mhoi'd^, 

which is precisely identity I5.4[ It is straightforward to check that P{f) is solution 
to l5.3[ One can easily show that T>{IV) is dense in Ker{r) which implies that: 



□ 



Vu G Ker{T), / {V P{f)Vu)dxdy = 0; 
Ju 

and therefore, one gets for all u in Ker{r) the following identity: 

/ \\/{P{f) + u)\'^dxdy= I \S/P{f)\^dxdy+ [ \\/u\^dxdy. 
Ju Ju Ju 

This entails that P{f) is the unique minimizer of 15.21 

Corollary 5.3. Let f be some element of H^^'^{M) and ]a,b[ an open interval . 
Consider the function m = f\]a,b[ and the associated variational problem: 

Suppose u is a solution of the following Dirichlet problem: 

'u G ^^(n) 

An = 

dn\9.\]a,b[ ^ 

Jiu\a,b[ = m 

Then T{u) is the unique solution to the above variational problem. 

Proof. Let g be an element of 7i^/^(]R) such that g\]a,b[ = Then we have the 
following identity: 
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I |2 



1 

1 



\V{P{g))\Hxdy 



V{u)\^dxdy 



{V{P{g) - u).Vu)dxdy 
> |r(n)|^i/2,^. + / {g{x) -T{u){x)) — {x,0)dx 

(Green) OU 



2n 



\V{P{g) - u)\^dxdy 



^ |r(M)|^i/2( 



5.2. Proof of theorem imi Let / be the trace on ] - 2L, 0[ of a 'h}/'^{M) function. 
By a general property of gaussian measures ([I3j), the conditional expectation on 
the left hand side can be obtained as the solution of the variational problem: 

By corollary 15.31 one must solve the Dirichlet problem 15.51 

'u e V^ill) 
Am = 
— =0 

an|R\]-2L,0[ 
.r(M)|]-2L,0[ = / 



The application u is solution to the above Dirichlet problem if and only if U given 
by U (x, y) = u{Lx — L, Ly) is solution to: 

(U eV^U) 
AU = 

an jR\]-l,l[ ^ 

[T{u)n-i,n = M 

where M{x) = f{Lx — L). One can solve this problem explicitly by using conformal 
mappings. Let (j) be the conformal map from the half disc 



{z e C; 1^1 < 1 andi?e(^) > 0} 



to the upper half plane H: 



1 



1, 



2i z' 

Then v{z) = U{(j){z)) is harmonic in D+, satisfies the condition |^ = on the vertical 
diameter {Re{z) = 0} fl and we have the following equality for x'^ + y'^ = 1: 



vix,y) = U{y,0) = M{y). 
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Since |^ = on the vertical diameter, one can extend symmetrically v as a harmonic 
function on the unit disc. By Poisson's formula, v has the following representation: 

^ ' 27r 1 - 2r cos(^ - t) + ^ ' 
In particular, one gets: 

1 W2 l_ 2 

v(0,y) = - -M(sin(t))dt. 

^ ' 7ry_,/2l-2ysin(t) + 2/2 ^ ^>> 

For X greater or equal to 1, 



t/(X,0) =v{<i)-\X,Q)) =t;(0,X- VX2-1). 
and for X less or equal to —1, 

U{X, 0) = t;(0, X + VX2 - 1). 

Therefore, one gets 

0) = - / ^ , , 7 --:r/(L(sin(t) - l))dt, 

where (/l is given by the formula: 



QLix) = I + 1 - + 1)2 - 1, if x > 0, 



or 



g^{x) = j + l + ^{^ + iy-l, ifx<-2L, 

One gets the desired result by making the change of variable s = L(sin(t) — 1) in 
the above integral. 

□ 

5.3. Proof of corollary 13.31 We take the limit L ^ oo in formula fl3.1|] . It is 

straightforward to check the convergence in 'D'QO, oo[) of the right hand term. Let 
us denote w1 the left hand term in (13.11) and w* the left hand term in (13. 2p . We 
have: 



and 



kllwi/^CR) = ,,,, inf kki/2( 



k*ki/2(R) = inf kki/2( 



Thus we have 7.^1/2 (k) ^ |tf *| 7^1/2 (jg). Modulo a subsequence, we may suppose 
by the Banach-Alaoglu theorem that w}^ converges weakly in n^/^{R) as L goes to 
infinity to w**. We have the convergence of w1 to w** in "D' \ which implies that 
^i]-oo,o[ = /• Since 7.^1/2 ^ | w*| 7^1/2 (k), we get w** = w*. 

5.4. Proof of theorem 13.51 
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The minimization problems. We consider the Sobolev space if^/^(M) equipped with 
the norm \.\t given by Ht'- 



Jr kr(t) 



where is the Fourier transform of kx'- 

kriO = ^sinc(27r|e|r). 

Let / be some function of ] — 2L,0[ which is the restriction of some function of 
ifi/2(R). We will denote by vl ^ the funct ion of iJ^/^(]R) solution of the minimization 
problem: 

inf Ivl^j. (Vt) 

We can consider the function / as a function of 7i^/^(]R) and note vl the unique 
solution of the minimization problem: 

A few intermediate propositions. We will start by giving a few intermediate propo- 
sitions before proving the theorem 13. 51 

Proposition 5.4. Let f be some function of H^^'^(K) .There exists a unique aL,T 
i7~^/^(M) with support in [— 2L, 0] that satisfies: 

Vt e] - 2L, 0[, (kr * «L,T)(t) = fit). (5.6) 

Proof. Consider the Hilbert space V of distributions in H~^^'^(K) with support in 
[— 2L,0]. On the product x V^, we consider the bilinear form: 



a{u,v) = I (kr * u){t)v{t)dt 
kriOmmdC 



By inequality (12. 2p . ^ a(u, u) defines a norm equivalent to iJ^^/^(]R). The linear 
form C : V Jj^ f{t)v{t)dt is continuous on V. By the Lax-Milgram theorem, there 
exists a unique aL,T ^ V such that: 

G V, a{aL,T,v) = C{v). 

Since Vi] — 2L, 0[) C V, we get the desired equality. □ 
In the sequel, we will denote by ipL^T the solution to equation (15.61) with / = 1 
and ip = (pi/2,1- Note that one can show the obvious scaling equality when 2L ^ T: 

^Mt)-- . (5.7) 

' 2L 1 + ln(T/2L) /_V(.sKs 
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with ip given by (cf. appendix): 

1 



VsG]-l,0[, v{s) 



27rln(2) ^-s-s^' 
We have thus the following expression for ^l,t'- 

1 1 



271 ln(2) + TT ln(T/2L) ^-s(s + 2L)' 



One can apply the above proposition to the resolution of the minimization prob- 
lems (Vt) and 

Proposition 5.5. Let ai.r be the solution of the equation ( (5. 6\) . Then the function 
v^T — ^T* (^L.T is solution to (Vt) ■ 

Proof Let v G H'^/^{R) be such that = / on ] - 2L, 0[. We get: 

Kt\t= I kT{i)\&LAi)?di 



v(t)aL,T(t)dt 
Similarly, wet get the following solution to (Voo): 



□ 



Proposition 5.6. Let f be in 7i-^/^(M). Let ai be the solution in H ^/^(M) with 
support in [— 2L, 0] of the equation: 

VtG]-2L,0[, (fc,^*a^)(t)=/(t)- ^-y/^ 'J)' J . 
Then v*t = k * + '^"^o ^5 solution to (Kxd). 

J-2L 'P(s/2L)ds 

Proof. First note that by the scaling relation (15.71) . we have: 

r_,Ms/2L)ds r_ 2LfL,2Lis)ds 
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Thus we get that: 



aL{t)dt= I {k2L*^L,2L){t)aL{t)dt 

{k2L * aL){t)'^L,2L{t)dt 



0. 



It is obvious by definition that ^2 is in 7Y^/^(M) and that f £ = / on the interval 
] — 2L, 0[. Now, if V satisfies the same conditions we get: 



Jr 

\M0\' 



m 



l I vl{t)aL{t)dt 



2 

_ 1 
~ 2 

^ K'|-Hl/2(IR)|w2|>^l/2 



v{t)aL{t)dt 



□ 

Remark 5.7. The solution is in fact not defined modulo constants and coincides 
with J_^^KL{t,s)f{s)ds for all f e H^/^{R) as can be seen by letting t go to or 
to infinity. 

Convergence of v*^ j. towards v*^ as T goes to infinity. We have the following exact 
formula: 

Proposition 5.8. LetT be larger than 2L. Then for allt in the interval [0,T — 2L], 
we have the identity: 



. , i-2Lm'Pi'mds i-(fc*y)(t/2L) 

' Ili^{s)ds 2L{l + \Yi{T/2L) j\^{s)ds) 
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Proof. For all t in the interval [0, T — 2L], we have since J^aL{s)ds = 0: 



fisMs/2L)ds 



vl{t)= / \n{T/\t-s\)aL{s)ds + 



2L /°2i<^(s/2L)ds 

/O pO 
ln+(r/|t-s|)K(s) -«L,T(s))rfs+ / ln^iT/\t-s\)aL,Tis)ds 



^ /_V fisMs/2L)ds 
f_^^^{s/2L)ds 



Now, we have for all t in the interval ] — 2L, 0[: 







f_,J{s)ip{s/2L)ds 



ln+(T/|t - s\){aL{s) - aL^T{s))ds 



-2L j_^^^[s/2L)ds 
By the uniqueness part of proposition 15.41 and the scahng relation (15. 7p . we get that: 



aiit) - aL,T{t) 



^{s/2L)ds 2L(1 + ln(T/2L) ^{s)ds) 



Plugging this relation into fl5.9l) yields the result. □ 
As a corollary of the above proposition, the explicit expression (13.31) and remark 
15. 7[ we easily get theorem 13.51 



6. Appendix 

In the appendix, we show that for all c ^ 0, the kernel ln'''(^;^) is positive (in 

particular, one can define a gaussian process with kernel (II. 3p ) and we show how to 
derive (13.31) . 

Lemma 6.1. The kernel hi^ {j^^ is positive. More precisely, we have the following 
expression for the Fourier transform if c < T: 
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Proof. We have for all ^: 
[ e-2^"«^ln+(^^)rfx = [ e-2^"«^ln(^^)cia; 

Jr fI + c J-(T-c) m + c 



(T-c) 

/O pT — C rji 

.(T-c) -X + C Jo X + C 

g— 2i7r5x pT—c ^~2i-K^x 

dx — — — - / dx 



J _{T_c) —X + c 2inC, Jq x + c 

1 f^~^ sin(2<a;) , 
dx 

■\ 

dx ^ 0. 



TTWo X + C 

1 sinfx 



7r|e|7o x + 27r|e|c 

□ 

In particular, one can define the gaussian process X'^ with covariance given by 
formula (11.31) . 

We now give a proof of identity (13. 3p . It is a direct consequence of the following 
lemma. 

Lemma 6.2. For allt in ]0,2[, we have: 

2 ds 
, ln|t-s|-== = -7rln(2) (6.1) 
V2s - 

Proof. We start by showing that the above quantity does not depend on t in the 
interval ]0,2[. Since the left hand side of (16.11) is invariant under t ^ 2 — t, we 
can consider the case t g]0, 1[. By performing successively the change of variables 
a = arcsin(l — s) and x = tan(cr/2), we get the following identities: 

2 



ds r''^ 

In |t — g| : = / In |t — 1 + sin(cr)|(icr 

V2s-s2 i„^/2 



1 r 

- / In |t — 1 + sin(o")|(i(T 

2 J — TV 

2x . dx 

In t - 1 



1 + x^ 1 + x^ 

1 + x'^ Jm 1 + x'^ 



For all e ^ 0, we introduce the following functions: 

Vte]0,2[, F,it) = l [ ln(|(l-t)(l + x2)-2xp + e)- 



dx 



x^ 



Note that we have to prove that Fq is a constant function on ]0, 1[. As e goes to 
zero, converges pointwise to Fq. We will prove that F^, tends to uniformly on 
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compact intervals [a, b] c]0, 1[ which imphes the result. If e > 0, we get: 

Let R be the rational fraction defined on C by: 

^ ' ((l-t)(l + z2)_2^)2 + e 

We get: 

1 (2;-2;+,+(t,0))(2-2_,_(t,0)) 



R{z) 



(1 - t) - e))(^ - e))(z - e))(^ - e)) ' 

where: 

By the formula of residues, we get that for all t in ]0, 1[: 

F'it) = ^ e) - z+,+ it, 0))(2:+,+ (t, e) - z.^.jt, 0)) 

1-t e) - z+^_{t, e))(2;+,+ (t, e) - e)){z+,+ {t, e) - e)) 

2m_ e) - z+,+ {t, 0))(^-,-(t, e) - 0)) 

1 - t e) - e))(2„,„(t, e) - z+^_{t, e))(2;_,_(t, e) - 2;„,+ (t, e)) ^" 

Using the above formula, one easily derives that F^' goes uniformly to on [a, b] as 
e goes to 0. 

Since Fq is a constant function on ]0, 2[, we only need to compute -Fo(l) to prove 
(16. ip . The computation of Fo{l) is standard. 

□ 
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Figure 1. First plot: empirical logvariogram of the SP500 index 
on the period 1989-1995. Second plot: empirical logvariogram of the 
SP500 index on the period 1995-2001. 
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Logvariogram SP500: 2001-2007 
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Logvariogram EUR/AUD: 2001-2007 
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Figure 2. First plot: empirical logvariogram of the SP500 index 
on the period 2001-2007. Second plot: empirical logvariogram of the 
currency EUR/AUD on the period 2001-2007 in the London time zone 
(Bloomberg series EUR/AUD CMPL) 
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